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ABSTRACT 

The cold dark matter (CDM) scenario for structure formation in the Universe 
is very attractive and has many successes; however, when its spectrum of 
density perturbations is normalized to the COBE anisotropy measurement 
the level of inhomogeneity predicted on small scales is too large. This can be 
remedied by a tau neutrino of mass 1 MeV — 10 MeV and lifetime 0.1 sec — 
100 sec whose decay products include electron neutrinos because it allows the 
total energy density in relativistic particles to be doubled without interfering 
with nucleosynthesis. The anisotropics predicted on the degree scale for 
"rCDM" are larger than standard CDM. Experiments at colliders may 
be able to probe such a mass range. 



The cold dark matter (CDM) scenario for the formation of structure in 
the Universe (i.e., galaxies, clusters of galaxies, superclusters, and so on) is 
both well motivated and very successful [0. The CDM model begins with a 
flat Universe with most of its mass in slowly moving ( "cold" ) particles such 
as axions or neutralinos and a baryon fraction compatible with primordial 
nucleosynthesis (~ 5%). The density perturbations that seed structure for- 
mation are the nearly scale-invariant perturbations that arise from quantum 
fluctuations during inflation . Standard CDM has one free parameter: the 
normalization of the spectrum. The recent COBE detection of anisotropy 
in the temperature of the cosmic background radiation (CBR) now provides 
the normalization by fixing the amplitude of perturbations on very large 
scales (A ~ lO^h'^ Mpc). But, when the spectrum is so normalized, the level 
of inhomogeneity predicted on small scales (A < 10h~^ Mpc) is too large by 
a factor of two (the Hubble constant Hq = 100/ikms~^ Mpc~^). Just how 
serious this problem is remains to be seen. 

A number of "fixes" have been proposed. They involve changing ei- 
ther the initial power spectrum or the energy content of the Universe. The 
simplest way of changing the power spectrum is "tilt" 0] (i.e., deviation 
from scale invariance) which reduces the power on small scales; in fact, 
several plausible models of inflation predict a tilted spectrum. However, 
this solution seems to lead to insufficient power on intermediate scales (A ~ 
30/i-^Mpc-lOO/i-^Mpc). Changing the energy content changes the transfer 
function, which relates the primeval power spectrum to that today p. For 
example, mixed dark matter (MDM) (65% CDM + 30% hot dark matter 
in the form of 7eV neutrinos + 5% baryons) or ACDM (80% vacuum 
energy -|- 15% CDM + 5% baryons) both reduce the power on small scales. 

The fix advocated here involves increasing the energy density in rela- 
tivistic particles, which also modifies the transfer function. In the stan- 
dard scenario the radiation content at late times {t ^ 1 sec) consists of 
photons and three massless neutrino species with slightly lower temper- 
ature = {A/liy^^Ty, accounting for a total radiation energy density 
Prad = g,TT^T^^/30. Here = 2 + 2(7/8) (4/11)^/3 jy^ counts the number 
of effectively massless degrees of freedom and is equal to 3.36 for = 3. In 
our "rCDM" model will be about a factor of two larger. 

To see why this helps, consider the scale Aeq, which crossed the horizon 
at the time the energy density in matter was equal to that in radiation. On 
scales much greater than Aeq the transfer function is unity because pertur- 
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bations on these scales crossed the horizon during the matter-dominated era 
and began growing as soon as they did. On scales much smaller than Aeq 
the transfer function decreases as because perturbations on these scales 
crossed the horizon before matter-radiation equality and did not begin grow- 
ing until the matter-dominated era. Thus, Aeq controls the shape of the 
transfer function; it is given by, Aeq ~ (-Ro/-Req)^eq {Rq is the scale factor 
today; Req and ^eq are the scale factor and age of the Universe at matter- 
radiation equality). Since the energy density in matter red shifts as i?"^ and 
that in radiation as Rq/Req = Pmatter/Prad — 2.4fio^^ x 10'^ /{g^/ 3. 36) 
(Qq is the ratio of the energy density to the critical density; f2o = 1 for a flat 
Universe). The age at matter-radiation equality ^eq ~ "^pi/ gl^'^T^Q, and so 

Aeq ~ 10 Mpc(no/i')"'(^7*/3.36)i/2_ 

Increasing g^, increases Aeq; since COBE fixes the power spectrum on scales 
A ^ Aeq, doing so reduces the power on small scales (see Fig. 1). 

How much additional radiation is needed? Since cosmological distances 
scale as the inverse of the Hubble constant, the shape of the power spectrum 
is set by (/i'^ Mpc)/Aeq oc T = {noh){gj3.3Q)-^/\ The best fit to all the 
data seems to require a "shape parameter" F ~ 0.3 whereas F ~ 0.5 
for the canonical values of Qq = 1-0, g* = 3.36, and h = 0.5^ For a given 
Hubble constant and shape parameter, the number of equivalent neutrino 
species called for is = 20{0.3/T)^nl{h/0.5y - 4.3. To achieve F = 0.3 
with h = 0.5(0.4) requires A^^ = 16(8.5). 

Therein lies the rub: Primordial nucleosynthesis restricts the equivalent 
number of massless neutrino species to be less than 3.4 [g^^ < 3.54) P]. The 
argument is well known: More relativistic degrees of freedom increase the 
expansion rate during nucleosynthesis and lead to an earlier freeze out of the 
neutron fraction at a higher value; this results in too much "^He if > 3.4. 

One way around this has been discussed previously: a late-decaying, mas- 
sive particle species, e.g., the once viable 17keV neutrino |T0|. Since the 
energy density in a massive particle species grows relative to radiation as the 
scale factor, a massive particle whose energy density at the epoch of nucle- 
osynthesis is negligible can produce a significant amount of radiation if it 



^In order to insure a sufficiently aged Universe with fio = 1, ^ must lay at the lower 
end of current measurements, h ~ 0.5. Of course, another fix is to posit a very low value 
of the Hubble parameter, h = 0.3, which leads to a very old Universe, to — 22 Gyr. 
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decays late enough, in the case of the 17keV neutrino, t ~ 10 ''sec |Tl]. 

Our proposal is very different and involves a subtle aspect of nucleosyn- 
thesis. We recently examined the effect of an unstable tau neutrino on nu- 
cleosynthesis for decay modes including 1/^ —* Uf. + (p {(f) is an effectively 
noninteracting scalar particle) [l^. The decay products of a massive tau 
neutrino can have more energy density than a massless neutrino species; for 
most decay modes this leads to an unacceptable increase in ^He. However, 
for the z/e0 decay mode there is a competing effect: Some decay-produced 
z/e's are captured by neutrons, reducing the neutron fraction and suppress- 
ing ^He production |]13|. This opens the door for adding relativistic degrees 
of freedom — which we've seen can help CDM — without overproducing "^He. 
The effect is even more pronounced for the i/^ — z/g + z/gt'e mode. 

The key to all of this is the fact that the decay-produced z/g's resemble 
a bath of neutrinos with temperature lower than that at which the neutron 
fraction freezes out (~ IMeV), and thus, their interactions with nucleons 
tend to reduce the neutron fraction. The rate for neutron-proton intercon- 
version (n + z/g ^ p + e~, etc.) varies as and so the typical energy of 
a captured neutrino is 5T, or 5MeV at freeze out; decay-produced electron 
neutrinos have energy (lower for the 3i>e mode). Thus, the effective 

temperature of decay-produced neutrinos will be lower than 1 MeV provided 
nil, ^ 10 MeV. Conversely, for masses greater than 10 MeV, decay-produced 
neutrinos increase the neutron fraction as they appear hotter. This simple 
explanation is borne out by detailed calculations [|l^ . 

Next, consider the lifetime dependence. Decays much earlier than 1 sec 
occur before freeze out, when the neutron fraction is close to its equilibrium 
value and so there is little effect. The dilution effect of the expansion de- 
creases the capture probability of a z/g, and so few z^e's produced much later 
than 1 sec are captured. Furthermore, the energy density of the massive tau 
neutrino grows relative to the radiation as the scale factor, and for lifetimes 
greater than 10 sec increases the expansion significantly rate, tending to in- 
crease '^He production. For long lifetimes the increase in expansion rate wins 
out and leads to higher "^He production; for lifetimes around a few seconds 
the two competing effects cancel; for short lifetimes the effect of neutron cap- 
tures wins out leading to lower ^He production. For a lifetime of ~ 0.3 sec 
the reduction in ^He is maximized and can be as large AYp ~ 0.10 [jT2| . 

There are three ways to take advantage of the "nucleosynthesis physics" 
described above to increase the radiation content in the Universe. The first 
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involves a tau neutrino of mass 1 MeV — 10 MeV and lifetime 0.1 sec — few sec. 
Tau-neutrino decays can reduce the neutron fraction enough to make room 
for up to the equivalent of 16 additional neutrino species (see Fig. 2). For 
lifetimes this short the energy density of the tau neutrino does not become 
significant before it decays, and the extra energy density in relativistic species 
must arise from additional massless degrees of freedom. Many extensions of 
the standard model predict new massless degrees of freedom (see below). 
(While we have not explored masses less than 1 MeV for technical reasons 
[|l^ , Fig. 2 suggests that a massive muon neutrino that decays to electron 
neutrinos might also allow for many additional massless species.) 

The second possibility involves a tau neutrino of lifetime 10 sec — 50 sec 
and is probably the most attractive. Here, not only do tau-neutrino decays 
inhibit ^He production but they also produce the additional energy density 
in relativistic particles. That energy density, expressed in equivalent number 
of neutrino species, is roughly 



AA/-^ ~ 0.4(rm^/MeV)(r^/sec)^/2^ (2) 

where r is the freeze-out abundance of the massive tau neutrino relative to 
that of a massless neutrino species (r ~ 0.1 — 1 for rrii, = 1 MeV — 10 MeV) 
and the total equivalent number of massless neutrinos Ni, = 2 + AN^. 

In Fig. 3 we show the regions of the tau-neutrino mass-lifetime plane that 
are excluded by primordial nucleosynthesis and contours of the energy den- 
sity of the tau-neutrino's decay products (based on our numerical calculations 
1^). We exclude a mass/lifetime if for no value of the baryon-to-photon ratio 



the light-element abundances satisfy: D/H > 10"^; (I)+^}ie)/B. < 1 x 10"^; 
Yp < 0.24; ''Li/H < 1.4 x 10^^°. As discussed earlier, this excludes long life- 
times. Since AA^,^ grows with lifetime, the boundary of the forbidden region 
is most interesting. Since it depends most sensitively on the limit to Yp (and 
to a lesser degree on ''Li) we also show the boundary for the relaxed con- 
straints, Yp < 0.25 and ^Li/H < 2 x 10~^°. For a Dirac neutrino AA^^, ~ 5(7) 
(i/e0) and 8(10) (Sz/g) are the maximum permitted (the number in parenthesis 
for the relaxed constraints). The AA"j,'s for a Majorana neutrino are about 



2 neutrino species smaller as their abundance is smaller fl^ . 

The final possibility is a variation on the previous one. While a massive 
tau neutrino provides a neat way of simultaneously suppressing "^He produc- 
tion and producing additional radiation, a species with similar decay modes 
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and a higher abundance could do even better. As noted earher, for a hfetime 
of a few seconds the effects of decay-produced z/g's and increased expansion 
rate cancel and increasing r does not affect ^He production — but does in- 
crease ANiy, cf. Eq. (0). Thus for r ~ 0.1 sec — few sec AA^,^ can be increased 
almost without limit. While r is not too different from one for a tau neutrino 
of mass 1 MeV — 10 MeV, it could be larger for a species with more degrees 
of freedom or weaker interactions. To illustrate, for a hypothetical 5 MeV 
species with r = 2, ANi, = 10 (i^e0), 17 (3z/e) as large as is permitted. 

Next, we consider a few of the consequences of rCDM. The level of 
inhomogeneity on small scales can be quantified by rms mass fluctuation 
in spheres of radius 8h~^ Mpc, ag- For COBE-normalized, standard CDM 
(jg — 1.2; for mixed dark matter ag — 0.75; and for rCDM model with 
r = 0.3, (jg ~ 0.67. The level of inhomogeneity predicted in rCDM on small 
scales is about a factor of two smaller than CDM and slightly lower than 
MDM. Since the rms fluctuation in galaxy counts is unity on the 8h~^ Mpc 
scale, light is a "biased" tracer of mass in both rCDM and MDM, with bias 
factor 6 = 1/(78 — 1-5. On the other hand, as can be seen in Fig. 1, rCDM 
has more power than MDM on very small scales and so forming objects such 
as quasars at red shifts z ~ 4 — 5 should be less problematic. 

We have computed the angular power spectrum of CBR anisotropy, Ci = 
(|a/mP) {ST = J^dimyim), and it is larger than standard CDM (see Fig. 4). 
To understand why, recall that the extra degrees of freedom delay matter- 
radiation equality. Last scattering occurs at the time of recombination 
(T ~ 1/3 eV), and in rCDM this happens closer to — or perhaps even in — 
the radiation-dominated era. As a result, the gravitational potential photons 
move through as they travel from the last-scattering surface is not constant. 



which leads to larger CBR anisotropy JTSl. The quantity l{l + l)Ci roughly 



corresponds to the expected variance of the CBR temperature difference mea- 
sured between directions separated by angle 9 ~ (200//)°. For I ~ 100 — 300 
this quantity is about a factor of 1.6 larger than standard CDM, which im- 
plies predicted temperature anisotropics on the degree scale that are a factor 
of a/1.6 larger. While the situation is far from settled, there are experiments 
with detections that are larger than the CDM prediction . 



In passing we mention two more consequences. First, a tau neutrino of 
mass 1 MeV — 10 MeV can resurrect a theory that predicts too many mass- 
less degrees of freedom to be consistent with the "standard nucleosynthesis 
bound," cf. Fig. 2. Second, the cosmic sea of i/g's is predicted to be very 
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different: a nonthermal spectrum with higher energy per neutrino. 

Finally, we turn to model building and verification. Within the framework 
of standard electroweak interactions the decay mode SUf, is highly forbidden, 
as it involves flavor-changing neutral currents, and of course the z/e0 mode 
does not exist. In simple extensions of the standard model both modes can 
arise. For example, in models with family symmetry is a massless Nambu- 
Goldstone boson associated with the breakdown of a family symmetry (at 
energy scale /) and ~ Sir f /ml ~ 10 sec(//109 GeV)2(m^/10 MeV)-^ 
In these models there can be many such Nambu-Goldstone bosons which 
contribute additional massless degress of freedom [0. The Sz/g decay mode 
can be mediated by a heavy boson with mass M; the lifetime is then t^, ~ 
0.1sec(m,y/10MeV)~^(M/Miy)^ sin~^ 9. We note though that in many models 
this boson would also couple to the tau lepton and violate the stringent 
bounds on r — > 3e. 

The best laboratory upper limits to the tau-neutrino mass are: 31 MeV 



by the ARGUS Collaboration and 32.6 MeV by the CLEO Collaboration p 



They are based upon end-point studies of the five-pion decay mode of the 
tau lepton. It is possible that the mass sensitivity could be extended to the 
5 MeV — 10 MeV range by looking at tau decay modes with multiple Kaons 



in the final state |T^; moreover, experiments done at B-factories may well 



be able to reach mass sensitivities of 5 MeV — 10 MeV |19 . 

To summarize, a tau neutrino of mass 1 MeV — 10 MeV and lifetime 
10 sec — 100 sec whose decay products include electron neutrinos can signifi- 
cantly increase the energy density in relativistic particles without upsetting 
the successful predictions of primordial nucleosynthesis. Alternatively, a tau 
neutrino in the same mass range, but with lifetime 0.1 sec — 1 sec, allows 
many more massless degrees of freedom leading to the same end. The higher 
level of energy density in relativistic particles modifies the CDM transfer 
function fixing the nagging problem that CDM has with excessive power on 
small scales. The level of CBR anisotropy predicted in the rCDM model is 
slightly larger than CDM on the degree scale. Most importantly, perhaps, 
prospects for testing this hypothesis at colliders appear promising. 
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Figure Captions 



Figure 1: The power spectra for CDM models normalized to COBE: stan- 
dard CDM (h = 0.5); mixed dark matter (h = 0.5); and rCDM {h = 0.45 and 
= 7.5). The data points are from the IRAS 1.2Jy survey |2y]. Wavenum- 
ber and wavelength are related hj k = 2tt/ X. 

Figure 2: Contours of the additional number of massless species that can be 
tolerated without violating the constraints to the light-element abundances 
(Dirac tau neutrino and Sz/g decay mode). [Note, because our calculations 
do not include inverse tau-neutrino decays, our results may be unreliable in 
the lower left-hand corner, r < ( MeV/m,^)^ sec; see Ref. |T2|.] 



Figure 3: The equivalent number of massless neutrino species AN^ pro- 
duced by the decay of a massive tau neutrino and the excluded regions of 
the mass-lifetime plane for a massive Dirac tau neutrino and the ^^(f) and 
3z/e modes. The solid contours correspond to AN^, = 4, 6, 8, 10 from left to 
right. The excluded regions are to the right of the broken lines (heavy for 
Yp < 0.24; light for Yp < 0.25). 

Figure 4: The COBE-normalized angular power spectrum for standard 
CDM (with h = 0.5) and rCDM (with h = .45 and AA^"^ = 10). For both, 
n^/i2 = 0.0125. 
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